\(\int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx\) [43]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (verified)
   Maple [A] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [A] (verification not implemented)
   Giac [B] (verification not implemented)
   Mupad [B] (verification not implemented)

Optimal result

Integrand size = 19, antiderivative size = 180 \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=\frac {4 a^3 b \text {arctanh}(\sin (c+d x))}{d}-\frac {6 a b^3 \text {arctanh}(\sin (c+d x))}{d}-\frac {a^4 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \cos (c+d x)}{d}-\frac {b^4 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \sec (c+d x)}{d}-\frac {2 b^4 \sec (c+d x)}{d}+\frac {b^4 \sec ^3(c+d x)}{3 d}-\frac {4 a^3 b \sin (c+d x)}{d}+\frac {6 a b^3 \sin (c+d x)}{d}+\frac {2 a b^3 \sin (c+d x) \tan ^2(c+d x)}{d} \]

[Out]

4*a^3*b*arctanh(sin(d*x+c))/d-6*a*b^3*arctanh(sin(d*x+c))/d-a^4*cos(d*x+c)/d+6*a^2*b^2*cos(d*x+c)/d-b^4*cos(d*
x+c)/d+6*a^2*b^2*sec(d*x+c)/d-2*b^4*sec(d*x+c)/d+1/3*b^4*sec(d*x+c)^3/d-4*a^3*b*sin(d*x+c)/d+6*a*b^3*sin(d*x+c
)/d+2*a*b^3*sin(d*x+c)*tan(d*x+c)^2/d

Rubi [A] (verified)

Time = 0.18 (sec) , antiderivative size = 180, normalized size of antiderivative = 1.00, number of steps used = 16, number of rules used = 9, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.474, Rules used = {3598, 2718, 2672, 327, 212, 2670, 14, 294, 276} \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=-\frac {a^4 \cos (c+d x)}{d}+\frac {4 a^3 b \text {arctanh}(\sin (c+d x))}{d}-\frac {4 a^3 b \sin (c+d x)}{d}+\frac {6 a^2 b^2 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \sec (c+d x)}{d}-\frac {6 a b^3 \text {arctanh}(\sin (c+d x))}{d}+\frac {6 a b^3 \sin (c+d x)}{d}+\frac {2 a b^3 \sin (c+d x) \tan ^2(c+d x)}{d}-\frac {b^4 \cos (c+d x)}{d}+\frac {b^4 \sec ^3(c+d x)}{3 d}-\frac {2 b^4 \sec (c+d x)}{d} \]

[In]

Int[Sin[c + d*x]*(a + b*Tan[c + d*x])^4,x]

[Out]

(4*a^3*b*ArcTanh[Sin[c + d*x]])/d - (6*a*b^3*ArcTanh[Sin[c + d*x]])/d - (a^4*Cos[c + d*x])/d + (6*a^2*b^2*Cos[
c + d*x])/d - (b^4*Cos[c + d*x])/d + (6*a^2*b^2*Sec[c + d*x])/d - (2*b^4*Sec[c + d*x])/d + (b^4*Sec[c + d*x]^3
)/(3*d) - (4*a^3*b*Sin[c + d*x])/d + (6*a*b^3*Sin[c + d*x])/d + (2*a*b^3*Sin[c + d*x]*Tan[c + d*x]^2)/d

Rule 14

Int[(u_)*((c_.)*(x_))^(m_.), x_Symbol] :> Int[ExpandIntegrand[(c*x)^m*u, x], x] /; FreeQ[{c, m}, x] && SumQ[u]
 &&  !LinearQ[u, x] &&  !MatchQ[u, (a_) + (b_.)*(v_) /; FreeQ[{a, b}, x] && InverseFunctionQ[v]]

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 276

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.), x_Symbol] :> Int[ExpandIntegrand[(c*x)^m*(a + b*x^n)^p,
 x], x] /; FreeQ[{a, b, c, m, n}, x] && IGtQ[p, 0]

Rule 294

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a + b*x^
n)^(p + 1)/(b*n*(p + 1))), x] - Dist[c^n*((m - n + 1)/(b*n*(p + 1))), Int[(c*x)^(m - n)*(a + b*x^n)^(p + 1), x
], x] /; FreeQ[{a, b, c}, x] && IGtQ[n, 0] && LtQ[p, -1] && GtQ[m + 1, n] &&  !ILtQ[(m + n*(p + 1) + 1)/n, 0]
&& IntBinomialQ[a, b, c, n, m, p, x]

Rule 327

Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Simp[c^(n - 1)*(c*x)^(m - n + 1)*((a + b*x^n
)^(p + 1)/(b*(m + n*p + 1))), x] - Dist[a*c^n*((m - n + 1)/(b*(m + n*p + 1))), Int[(c*x)^(m - n)*(a + b*x^n)^p
, x], x] /; FreeQ[{a, b, c, p}, x] && IGtQ[n, 0] && GtQ[m, n - 1] && NeQ[m + n*p + 1, 0] && IntBinomialQ[a, b,
 c, n, m, p, x]

Rule 2670

Int[sin[(e_.) + (f_.)*(x_)]^(m_.)*tan[(e_.) + (f_.)*(x_)]^(n_.), x_Symbol] :> Dist[-f^(-1), Subst[Int[(1 - x^2
)^((m + n - 1)/2)/x^n, x], x, Cos[e + f*x]], x] /; FreeQ[{e, f}, x] && IntegersQ[m, n, (m + n - 1)/2]

Rule 2672

Int[((a_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*tan[(e_.) + (f_.)*(x_)]^(n_.), x_Symbol] :> With[{ff = FreeFactors[S
in[e + f*x], x]}, Dist[ff/f, Subst[Int[(ff*x)^(m + n)/(a^2 - ff^2*x^2)^((n + 1)/2), x], x, a*(Sin[e + f*x]/ff)
], x]] /; FreeQ[{a, e, f, m}, x] && IntegerQ[(n + 1)/2]

Rule 2718

Int[sin[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-Cos[c + d*x]/d, x] /; FreeQ[{c, d}, x]

Rule 3598

Int[sin[(e_.) + (f_.)*(x_)]^(m_.)*((a_) + (b_.)*tan[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Int[Expand[Sin[e
+ f*x]^m*(a + b*Tan[e + f*x])^n, x], x] /; FreeQ[{a, b, e, f}, x] && IntegerQ[(m - 1)/2] && IGtQ[n, 0]

Rubi steps \begin{align*} \text {integral}& = \int \left (a^4 \sin (c+d x)+4 a^3 b \sin (c+d x) \tan (c+d x)+6 a^2 b^2 \sin (c+d x) \tan ^2(c+d x)+4 a b^3 \sin (c+d x) \tan ^3(c+d x)+b^4 \sin (c+d x) \tan ^4(c+d x)\right ) \, dx \\ & = a^4 \int \sin (c+d x) \, dx+\left (4 a^3 b\right ) \int \sin (c+d x) \tan (c+d x) \, dx+\left (6 a^2 b^2\right ) \int \sin (c+d x) \tan ^2(c+d x) \, dx+\left (4 a b^3\right ) \int \sin (c+d x) \tan ^3(c+d x) \, dx+b^4 \int \sin (c+d x) \tan ^4(c+d x) \, dx \\ & = -\frac {a^4 \cos (c+d x)}{d}+\frac {\left (4 a^3 b\right ) \text {Subst}\left (\int \frac {x^2}{1-x^2} \, dx,x,\sin (c+d x)\right )}{d}-\frac {\left (6 a^2 b^2\right ) \text {Subst}\left (\int \frac {1-x^2}{x^2} \, dx,x,\cos (c+d x)\right )}{d}+\frac {\left (4 a b^3\right ) \text {Subst}\left (\int \frac {x^4}{\left (1-x^2\right )^2} \, dx,x,\sin (c+d x)\right )}{d}-\frac {b^4 \text {Subst}\left (\int \frac {\left (1-x^2\right )^2}{x^4} \, dx,x,\cos (c+d x)\right )}{d} \\ & = -\frac {a^4 \cos (c+d x)}{d}-\frac {4 a^3 b \sin (c+d x)}{d}+\frac {2 a b^3 \sin (c+d x) \tan ^2(c+d x)}{d}+\frac {\left (4 a^3 b\right ) \text {Subst}\left (\int \frac {1}{1-x^2} \, dx,x,\sin (c+d x)\right )}{d}-\frac {\left (6 a^2 b^2\right ) \text {Subst}\left (\int \left (-1+\frac {1}{x^2}\right ) \, dx,x,\cos (c+d x)\right )}{d}-\frac {\left (6 a b^3\right ) \text {Subst}\left (\int \frac {x^2}{1-x^2} \, dx,x,\sin (c+d x)\right )}{d}-\frac {b^4 \text {Subst}\left (\int \left (1+\frac {1}{x^4}-\frac {2}{x^2}\right ) \, dx,x,\cos (c+d x)\right )}{d} \\ & = \frac {4 a^3 b \text {arctanh}(\sin (c+d x))}{d}-\frac {a^4 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \cos (c+d x)}{d}-\frac {b^4 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \sec (c+d x)}{d}-\frac {2 b^4 \sec (c+d x)}{d}+\frac {b^4 \sec ^3(c+d x)}{3 d}-\frac {4 a^3 b \sin (c+d x)}{d}+\frac {6 a b^3 \sin (c+d x)}{d}+\frac {2 a b^3 \sin (c+d x) \tan ^2(c+d x)}{d}-\frac {\left (6 a b^3\right ) \text {Subst}\left (\int \frac {1}{1-x^2} \, dx,x,\sin (c+d x)\right )}{d} \\ & = \frac {4 a^3 b \text {arctanh}(\sin (c+d x))}{d}-\frac {6 a b^3 \text {arctanh}(\sin (c+d x))}{d}-\frac {a^4 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \cos (c+d x)}{d}-\frac {b^4 \cos (c+d x)}{d}+\frac {6 a^2 b^2 \sec (c+d x)}{d}-\frac {2 b^4 \sec (c+d x)}{d}+\frac {b^4 \sec ^3(c+d x)}{3 d}-\frac {4 a^3 b \sin (c+d x)}{d}+\frac {6 a b^3 \sin (c+d x)}{d}+\frac {2 a b^3 \sin (c+d x) \tan ^2(c+d x)}{d} \\ \end{align*}

Mathematica [B] (verified)

Leaf count is larger than twice the leaf count of optimal. \(383\) vs. \(2(180)=360\).

Time = 6.41 (sec) , antiderivative size = 383, normalized size of antiderivative = 2.13 \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=\frac {72 a^2 b^2-22 b^4-12 \left (a^4-6 a^2 b^2+b^4\right ) \cos (c+d x)-24 a b \left (2 a^2-3 b^2\right ) \log \left (\cos \left (\frac {1}{2} (c+d x)\right )-\sin \left (\frac {1}{2} (c+d x)\right )\right )+24 a b \left (2 a^2-3 b^2\right ) \log \left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right )+\frac {b^3 (12 a+b)}{\left (\cos \left (\frac {1}{2} (c+d x)\right )-\sin \left (\frac {1}{2} (c+d x)\right )\right )^2}+\frac {2 b^4 \sin \left (\frac {1}{2} (c+d x)\right )}{\left (\cos \left (\frac {1}{2} (c+d x)\right )-\sin \left (\frac {1}{2} (c+d x)\right )\right )^3}+\frac {2 b^2 \left (36 a^2-11 b^2\right ) \sin \left (\frac {1}{2} (c+d x)\right )}{\cos \left (\frac {1}{2} (c+d x)\right )-\sin \left (\frac {1}{2} (c+d x)\right )}-\frac {2 b^4 \sin \left (\frac {1}{2} (c+d x)\right )}{\left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right )^3}+\frac {b^3 (-12 a+b)}{\left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right )^2}+\frac {2 b^2 \left (-36 a^2+11 b^2\right ) \sin \left (\frac {1}{2} (c+d x)\right )}{\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )}-48 a b \left (a^2-b^2\right ) \sin (c+d x)}{12 d} \]

[In]

Integrate[Sin[c + d*x]*(a + b*Tan[c + d*x])^4,x]

[Out]

(72*a^2*b^2 - 22*b^4 - 12*(a^4 - 6*a^2*b^2 + b^4)*Cos[c + d*x] - 24*a*b*(2*a^2 - 3*b^2)*Log[Cos[(c + d*x)/2] -
 Sin[(c + d*x)/2]] + 24*a*b*(2*a^2 - 3*b^2)*Log[Cos[(c + d*x)/2] + Sin[(c + d*x)/2]] + (b^3*(12*a + b))/(Cos[(
c + d*x)/2] - Sin[(c + d*x)/2])^2 + (2*b^4*Sin[(c + d*x)/2])/(Cos[(c + d*x)/2] - Sin[(c + d*x)/2])^3 + (2*b^2*
(36*a^2 - 11*b^2)*Sin[(c + d*x)/2])/(Cos[(c + d*x)/2] - Sin[(c + d*x)/2]) - (2*b^4*Sin[(c + d*x)/2])/(Cos[(c +
 d*x)/2] + Sin[(c + d*x)/2])^3 + (b^3*(-12*a + b))/(Cos[(c + d*x)/2] + Sin[(c + d*x)/2])^2 + (2*b^2*(-36*a^2 +
 11*b^2)*Sin[(c + d*x)/2])/(Cos[(c + d*x)/2] + Sin[(c + d*x)/2]) - 48*a*b*(a^2 - b^2)*Sin[c + d*x])/(12*d)

Maple [A] (verified)

Time = 3.04 (sec) , antiderivative size = 217, normalized size of antiderivative = 1.21

method result size
derivativedivides \(\frac {-a^{4} \cos \left (d x +c \right )+4 a^{3} b \left (-\sin \left (d x +c \right )+\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )\right )+6 a^{2} b^{2} \left (\frac {\sin ^{4}\left (d x +c \right )}{\cos \left (d x +c \right )}+\left (2+\sin ^{2}\left (d x +c \right )\right ) \cos \left (d x +c \right )\right )+4 a \,b^{3} \left (\frac {\sin ^{5}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{3}\left (d x +c \right )\right )}{2}+\frac {3 \sin \left (d x +c \right )}{2}-\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+b^{4} \left (\frac {\sin ^{6}\left (d x +c \right )}{3 \cos \left (d x +c \right )^{3}}-\frac {\sin ^{6}\left (d x +c \right )}{\cos \left (d x +c \right )}-\left (\frac {8}{3}+\sin ^{4}\left (d x +c \right )+\frac {4 \left (\sin ^{2}\left (d x +c \right )\right )}{3}\right ) \cos \left (d x +c \right )\right )}{d}\) \(217\)
default \(\frac {-a^{4} \cos \left (d x +c \right )+4 a^{3} b \left (-\sin \left (d x +c \right )+\ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )\right )+6 a^{2} b^{2} \left (\frac {\sin ^{4}\left (d x +c \right )}{\cos \left (d x +c \right )}+\left (2+\sin ^{2}\left (d x +c \right )\right ) \cos \left (d x +c \right )\right )+4 a \,b^{3} \left (\frac {\sin ^{5}\left (d x +c \right )}{2 \cos \left (d x +c \right )^{2}}+\frac {\left (\sin ^{3}\left (d x +c \right )\right )}{2}+\frac {3 \sin \left (d x +c \right )}{2}-\frac {3 \ln \left (\sec \left (d x +c \right )+\tan \left (d x +c \right )\right )}{2}\right )+b^{4} \left (\frac {\sin ^{6}\left (d x +c \right )}{3 \cos \left (d x +c \right )^{3}}-\frac {\sin ^{6}\left (d x +c \right )}{\cos \left (d x +c \right )}-\left (\frac {8}{3}+\sin ^{4}\left (d x +c \right )+\frac {4 \left (\sin ^{2}\left (d x +c \right )\right )}{3}\right ) \cos \left (d x +c \right )\right )}{d}\) \(217\)
risch \(\frac {2 i {\mathrm e}^{i \left (d x +c \right )} a^{3} b}{d}-\frac {2 i {\mathrm e}^{i \left (d x +c \right )} a \,b^{3}}{d}-\frac {{\mathrm e}^{i \left (d x +c \right )} a^{4}}{2 d}+\frac {3 \,{\mathrm e}^{i \left (d x +c \right )} a^{2} b^{2}}{d}-\frac {{\mathrm e}^{i \left (d x +c \right )} b^{4}}{2 d}-\frac {2 i {\mathrm e}^{-i \left (d x +c \right )} a^{3} b}{d}+\frac {2 i {\mathrm e}^{-i \left (d x +c \right )} a \,b^{3}}{d}-\frac {{\mathrm e}^{-i \left (d x +c \right )} a^{4}}{2 d}+\frac {3 \,{\mathrm e}^{-i \left (d x +c \right )} a^{2} b^{2}}{d}-\frac {{\mathrm e}^{-i \left (d x +c \right )} b^{4}}{2 d}-\frac {4 b^{2} {\mathrm e}^{i \left (d x +c \right )} \left (-9 a^{2} {\mathrm e}^{4 i \left (d x +c \right )}+3 b^{2} {\mathrm e}^{4 i \left (d x +c \right )}+3 i a b \,{\mathrm e}^{4 i \left (d x +c \right )}-18 a^{2} {\mathrm e}^{2 i \left (d x +c \right )}+4 b^{2} {\mathrm e}^{2 i \left (d x +c \right )}-9 a^{2}+3 b^{2}-3 i a b \right )}{3 d \left ({\mathrm e}^{2 i \left (d x +c \right )}+1\right )^{3}}+\frac {4 a^{3} b \ln \left ({\mathrm e}^{i \left (d x +c \right )}+i\right )}{d}-\frac {6 a \,b^{3} \ln \left ({\mathrm e}^{i \left (d x +c \right )}+i\right )}{d}-\frac {4 a^{3} b \ln \left ({\mathrm e}^{i \left (d x +c \right )}-i\right )}{d}+\frac {6 a \,b^{3} \ln \left ({\mathrm e}^{i \left (d x +c \right )}-i\right )}{d}\) \(412\)

[In]

int(sin(d*x+c)*(a+b*tan(d*x+c))^4,x,method=_RETURNVERBOSE)

[Out]

1/d*(-a^4*cos(d*x+c)+4*a^3*b*(-sin(d*x+c)+ln(sec(d*x+c)+tan(d*x+c)))+6*a^2*b^2*(sin(d*x+c)^4/cos(d*x+c)+(2+sin
(d*x+c)^2)*cos(d*x+c))+4*a*b^3*(1/2*sin(d*x+c)^5/cos(d*x+c)^2+1/2*sin(d*x+c)^3+3/2*sin(d*x+c)-3/2*ln(sec(d*x+c
)+tan(d*x+c)))+b^4*(1/3*sin(d*x+c)^6/cos(d*x+c)^3-sin(d*x+c)^6/cos(d*x+c)-(8/3+sin(d*x+c)^4+4/3*sin(d*x+c)^2)*
cos(d*x+c)))

Fricas [A] (verification not implemented)

none

Time = 0.30 (sec) , antiderivative size = 176, normalized size of antiderivative = 0.98 \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=-\frac {3 \, {\left (a^{4} - 6 \, a^{2} b^{2} + b^{4}\right )} \cos \left (d x + c\right )^{4} - 3 \, {\left (2 \, a^{3} b - 3 \, a b^{3}\right )} \cos \left (d x + c\right )^{3} \log \left (\sin \left (d x + c\right ) + 1\right ) + 3 \, {\left (2 \, a^{3} b - 3 \, a b^{3}\right )} \cos \left (d x + c\right )^{3} \log \left (-\sin \left (d x + c\right ) + 1\right ) - b^{4} - 6 \, {\left (3 \, a^{2} b^{2} - b^{4}\right )} \cos \left (d x + c\right )^{2} - 6 \, {\left (a b^{3} \cos \left (d x + c\right ) - 2 \, {\left (a^{3} b - a b^{3}\right )} \cos \left (d x + c\right )^{3}\right )} \sin \left (d x + c\right )}{3 \, d \cos \left (d x + c\right )^{3}} \]

[In]

integrate(sin(d*x+c)*(a+b*tan(d*x+c))^4,x, algorithm="fricas")

[Out]

-1/3*(3*(a^4 - 6*a^2*b^2 + b^4)*cos(d*x + c)^4 - 3*(2*a^3*b - 3*a*b^3)*cos(d*x + c)^3*log(sin(d*x + c) + 1) +
3*(2*a^3*b - 3*a*b^3)*cos(d*x + c)^3*log(-sin(d*x + c) + 1) - b^4 - 6*(3*a^2*b^2 - b^4)*cos(d*x + c)^2 - 6*(a*
b^3*cos(d*x + c) - 2*(a^3*b - a*b^3)*cos(d*x + c)^3)*sin(d*x + c))/(d*cos(d*x + c)^3)

Sympy [F]

\[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=\int \left (a + b \tan {\left (c + d x \right )}\right )^{4} \sin {\left (c + d x \right )}\, dx \]

[In]

integrate(sin(d*x+c)*(a+b*tan(d*x+c))**4,x)

[Out]

Integral((a + b*tan(c + d*x))**4*sin(c + d*x), x)

Maxima [A] (verification not implemented)

none

Time = 0.21 (sec) , antiderivative size = 166, normalized size of antiderivative = 0.92 \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=-\frac {3 \, a b^{3} {\left (\frac {2 \, \sin \left (d x + c\right )}{\sin \left (d x + c\right )^{2} - 1} + 3 \, \log \left (\sin \left (d x + c\right ) + 1\right ) - 3 \, \log \left (\sin \left (d x + c\right ) - 1\right ) - 4 \, \sin \left (d x + c\right )\right )} - 18 \, a^{2} b^{2} {\left (\frac {1}{\cos \left (d x + c\right )} + \cos \left (d x + c\right )\right )} + b^{4} {\left (\frac {6 \, \cos \left (d x + c\right )^{2} - 1}{\cos \left (d x + c\right )^{3}} + 3 \, \cos \left (d x + c\right )\right )} - 6 \, a^{3} b {\left (\log \left (\sin \left (d x + c\right ) + 1\right ) - \log \left (\sin \left (d x + c\right ) - 1\right ) - 2 \, \sin \left (d x + c\right )\right )} + 3 \, a^{4} \cos \left (d x + c\right )}{3 \, d} \]

[In]

integrate(sin(d*x+c)*(a+b*tan(d*x+c))^4,x, algorithm="maxima")

[Out]

-1/3*(3*a*b^3*(2*sin(d*x + c)/(sin(d*x + c)^2 - 1) + 3*log(sin(d*x + c) + 1) - 3*log(sin(d*x + c) - 1) - 4*sin
(d*x + c)) - 18*a^2*b^2*(1/cos(d*x + c) + cos(d*x + c)) + b^4*((6*cos(d*x + c)^2 - 1)/cos(d*x + c)^3 + 3*cos(d
*x + c)) - 6*a^3*b*(log(sin(d*x + c) + 1) - log(sin(d*x + c) - 1) - 2*sin(d*x + c)) + 3*a^4*cos(d*x + c))/d

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 19074 vs. \(2 (178) = 356\).

Time = 11.22 (sec) , antiderivative size = 19074, normalized size of antiderivative = 105.97 \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=\text {Too large to display} \]

[In]

integrate(sin(d*x+c)*(a+b*tan(d*x+c))^4,x, algorithm="giac")

[Out]

-1/3*(6*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 +
 tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x
)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^8 - 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*
d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c
) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^8 - 6*a^3*
b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x
)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/
2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^8 + 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1
/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(
1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^8 + 3*a^4*tan(1/2*d*x)
^8*tan(1/2*c)^8 - 36*a^2*b^2*tan(1/2*d*x)^8*tan(1/2*c)^8 + 8*b^4*tan(1/2*d*x)^8*tan(1/2*c)^8 - 12*a^3*b*log(2*
(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + ta
n(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 +
 1))*tan(1/2*d*x)^8*tan(1/2*c)^6 + 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) +
 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x
)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^6 + 12*a^3*b*log(2*(tan(1/2*d
*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2
 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1
/2*d*x)^8*tan(1/2*c)^6 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2
*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/
2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^6 - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(
1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1
/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*
tan(1/2*c)^7 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan
(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 +
tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^7 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2
- 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) +
 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c
)^7 - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2
 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d
*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^7 - 24*a^3*b*tan(1/2*d*x)^8*tan(1/2*c)^7 + 36*a*b^3*tan(1
/2*d*x)^8*tan(1/2*c)^7 - 12*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2
*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/
2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^8 + 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(
1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1
/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*
tan(1/2*c)^8 + 12*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(
1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + t
an(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^8 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 -
 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) +
2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)
^8 - 24*a^3*b*tan(1/2*d*x)^7*tan(1/2*c)^8 + 36*a*b^3*tan(1/2*d*x)^7*tan(1/2*c)^8 - 12*a^4*tan(1/2*d*x)^8*tan(1
/2*c)^6 + 72*a^2*b^2*tan(1/2*d*x)^8*tan(1/2*c)^6 - 16*b^4*tan(1/2*d*x)^8*tan(1/2*c)^6 - 48*a^4*tan(1/2*d*x)^7*
tan(1/2*c)^7 + 432*a^2*b^2*tan(1/2*d*x)^7*tan(1/2*c)^7 - 96*b^4*tan(1/2*d*x)^7*tan(1/2*c)^7 - 12*a^4*tan(1/2*d
*x)^6*tan(1/2*c)^8 + 72*a^2*b^2*tan(1/2*d*x)^6*tan(1/2*c)^8 - 16*b^4*tan(1/2*d*x)^6*tan(1/2*c)^8 + 72*a^3*b*lo
g(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2
+ tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)
^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^5 - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2
*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/
2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^5 - 72*a^3*b*log(2*(tan(
1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*
tan(1/2*d*x)^7*tan(1/2*c)^5 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*t
an(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*
tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^5 + 72*a^3*b*tan(1/2*d*x)^8*tan(1
/2*c)^5 - 60*a*b^3*tan(1/2*d*x)^8*tan(1/2*c)^5 + 312*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)
^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) +
1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^6 - 468*a*b^3*
log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^
2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*
c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^6 - 312*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1
/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(
1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^6 + 468*a*b^3*log(2*(t
an(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(
1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1
))*tan(1/2*d*x)^6*tan(1/2*c)^6 + 384*a^3*b*tan(1/2*d*x)^7*tan(1/2*c)^6 - 432*a*b^3*tan(1/2*d*x)^7*tan(1/2*c)^6
 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + t
an(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^
2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^7 - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*
d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c
) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^7 - 72*a^3
*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*
x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1
/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^7 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*ta
n(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(t
an(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^7 + 384*a^3*b*tan(1
/2*d*x)^6*tan(1/2*c)^7 - 432*a*b^3*tan(1/2*d*x)^6*tan(1/2*c)^7 + 72*a^3*b*tan(1/2*d*x)^5*tan(1/2*c)^8 - 60*a*b
^3*tan(1/2*d*x)^5*tan(1/2*c)^8 + 18*a^4*tan(1/2*d*x)^8*tan(1/2*c)^4 - 72*a^2*b^2*tan(1/2*d*x)^8*tan(1/2*c)^4 +
 144*a^4*tan(1/2*d*x)^7*tan(1/2*c)^5 - 720*a^2*b^2*tan(1/2*d*x)^7*tan(1/2*c)^5 + 96*b^4*tan(1/2*d*x)^7*tan(1/2
*c)^5 + 336*a^4*tan(1/2*d*x)^6*tan(1/2*c)^6 - 2304*a^2*b^2*tan(1/2*d*x)^6*tan(1/2*c)^6 + 416*b^4*tan(1/2*d*x)^
6*tan(1/2*c)^6 + 144*a^4*tan(1/2*d*x)^5*tan(1/2*c)^7 - 720*a^2*b^2*tan(1/2*d*x)^5*tan(1/2*c)^7 + 96*b^4*tan(1/
2*d*x)^5*tan(1/2*c)^7 + 18*a^4*tan(1/2*d*x)^4*tan(1/2*c)^8 - 72*a^2*b^2*tan(1/2*d*x)^4*tan(1/2*c)^8 + 12*a^3*b
*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^2 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1
/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(
1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^2 - 12*a^3*b*log(2*(ta
n(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1
/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1)
)*tan(1/2*d*x)^8*tan(1/2*c)^2 + 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*
tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8*tan(1/2*c)^2 + 72*a^3*b*log(2*(tan(1/2*d*x)
^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 -
2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*
d*x)^7*tan(1/2*c)^3 - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d
*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*
c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^3 - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/
2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2
*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*ta
n(1/2*c)^3 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1
/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + ta
n(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c)^3 - 72*a^3*b*tan(1/2*d*x)^8*tan(1/2*c)^3 + 60*a*b^
3*tan(1/2*d*x)^8*tan(1/2*c)^3 - 720*a^3*b*tan(1/2*d*x)^7*tan(1/2*c)^4 + 600*a*b^3*tan(1/2*d*x)^7*tan(1/2*c)^4
- 456*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + t
an(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^
2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^5 + 684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*
d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c
) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^5 + 456*a^
3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d
*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(
1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^5 - 684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*t
an(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(
tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^5 - 2304*a^3*b*tan
(1/2*d*x)^6*tan(1/2*c)^5 + 2256*a*b^3*tan(1/2*d*x)^6*tan(1/2*c)^5 - 2304*a^3*b*tan(1/2*d*x)^5*tan(1/2*c)^6 + 2
256*a*b^3*tan(1/2*d*x)^5*tan(1/2*c)^6 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2
*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/
2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^7 - 108*a*b^3*log(2*(tan
(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/
2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))
*tan(1/2*d*x)^3*tan(1/2*c)^7 - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*t
an(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*
tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^7 + 108*a*b^3*log(2*(tan(1/2*d*x)
^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 +
2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*
d*x)^3*tan(1/2*c)^7 - 720*a^3*b*tan(1/2*d*x)^4*tan(1/2*c)^7 + 600*a*b^3*tan(1/2*d*x)^4*tan(1/2*c)^7 + 12*a^3*b
*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^8 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1
/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(
1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^8 - 12*a^3*b*log(2*(ta
n(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1
/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1)
)*tan(1/2*d*x)^2*tan(1/2*c)^8 + 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*
tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^8 - 72*a^3*b*tan(1/2*d*x)^3*tan(
1/2*c)^8 + 60*a*b^3*tan(1/2*d*x)^3*tan(1/2*c)^8 - 12*a^4*tan(1/2*d*x)^8*tan(1/2*c)^2 + 72*a^2*b^2*tan(1/2*d*x)
^8*tan(1/2*c)^2 - 16*b^4*tan(1/2*d*x)^8*tan(1/2*c)^2 - 144*a^4*tan(1/2*d*x)^7*tan(1/2*c)^3 + 720*a^2*b^2*tan(1
/2*d*x)^7*tan(1/2*c)^3 - 96*b^4*tan(1/2*d*x)^7*tan(1/2*c)^3 - 648*a^4*tan(1/2*d*x)^6*tan(1/2*c)^4 + 3312*a^2*b
^2*tan(1/2*d*x)^6*tan(1/2*c)^4 - 480*b^4*tan(1/2*d*x)^6*tan(1/2*c)^4 - 1200*a^4*tan(1/2*d*x)^5*tan(1/2*c)^5 +
7344*a^2*b^2*tan(1/2*d*x)^5*tan(1/2*c)^5 - 1248*b^4*tan(1/2*d*x)^5*tan(1/2*c)^5 - 648*a^4*tan(1/2*d*x)^4*tan(1
/2*c)^6 + 3312*a^2*b^2*tan(1/2*d*x)^4*tan(1/2*c)^6 - 480*b^4*tan(1/2*d*x)^4*tan(1/2*c)^6 - 144*a^4*tan(1/2*d*x
)^3*tan(1/2*c)^7 + 720*a^2*b^2*tan(1/2*d*x)^3*tan(1/2*c)^7 - 96*b^4*tan(1/2*d*x)^3*tan(1/2*c)^7 - 12*a^4*tan(1
/2*d*x)^2*tan(1/2*c)^8 + 72*a^2*b^2*tan(1/2*d*x)^2*tan(1/2*c)^8 - 16*b^4*tan(1/2*d*x)^2*tan(1/2*c)^8 - 6*a^3*b
*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 1))*tan(1/2*d*x)^8 + 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(
1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan
(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8 + 6*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 -
2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2
*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^8 - 9*a*b^3*l
og(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2
 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c
)^2 + 1))*tan(1/2*d*x)^8 - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1
/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(
1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c) + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan
(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(
1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7
*tan(1/2*c) + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1
/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + ta
n(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c) - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2
*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*
tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^7*tan(1/2*c) +
 24*a^3*b*tan(1/2*d*x)^8*tan(1/2*c) - 36*a*b^3*tan(1/2*d*x)^8*tan(1/2*c) - 312*a^3*b*log(2*(tan(1/2*d*x)^2*tan
(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(
1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6
*tan(1/2*c)^2 + 468*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*ta
n(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 +
 tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^2 + 312*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^
2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x)
 + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2
*c)^2 - 468*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6*tan(1/2*c)^2 + 384*a^3*b*tan(1/2*d*x)^7*tan(1/2*c)^2 - 432*a*b^3*t
an(1/2*d*x)^7*tan(1/2*c)^2 - 456*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*ta
n(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*t
an(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^3 + 684*a*b^3*log(2*(tan(1/2*d*x)^
2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2
*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d
*x)^5*tan(1/2*c)^3 + 456*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*
x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c
)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c)^3 - 684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/
2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2
*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*ta
n(1/2*c)^3 + 2304*a^3*b*tan(1/2*d*x)^6*tan(1/2*c)^3 - 2256*a*b^3*tan(1/2*d*x)^6*tan(1/2*c)^3 + 6000*a^3*b*tan(
1/2*d*x)^5*tan(1/2*c)^4 - 6120*a*b^3*tan(1/2*d*x)^5*tan(1/2*c)^4 - 456*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^
2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x)
 - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2
*c)^5 + 684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^5 + 456*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*ta
n(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan
(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^5 -
684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan
(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2
+ tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^5 + 6000*a^3*b*tan(1/2*d*x)^4*tan(1/2*c)^5 - 6120*a*b^3*tan(1/2
*d*x)^4*tan(1/2*c)^5 - 312*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*
d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2
*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^6 + 468*a*b^3*log(2*(tan(1/2*d*x)^2*tan(
1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1
/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*
tan(1/2*c)^6 + 312*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan
(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 +
tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^6 - 468*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2
 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x)
+ 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*
c)^6 + 2304*a^3*b*tan(1/2*d*x)^3*tan(1/2*c)^6 - 2256*a*b^3*tan(1/2*d*x)^3*tan(1/2*c)^6 - 72*a^3*b*log(2*(tan(1
/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*
c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*t
an(1/2*d*x)*tan(1/2*c)^7 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(
1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan
(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^7 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan
(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(
1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*t
an(1/2*c)^7 - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(
1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + t
an(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^7 + 384*a^3*b*tan(1/2*d*x)^2*tan(1/2*c)^7 - 432*a*b
^3*tan(1/2*d*x)^2*tan(1/2*c)^7 - 6*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*
tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^8 + 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2
+ 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) -
 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^8 + 6*a^3*b*l
og(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2
 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c
)^2 + 1))*tan(1/2*c)^8 - 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*
d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2
*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^8 + 24*a^3*b*tan(1/2*d*x)*tan(1/2*c)^8 - 36*a*b^3*tan(1
/2*d*x)*tan(1/2*c)^8 + 3*a^4*tan(1/2*d*x)^8 - 36*a^2*b^2*tan(1/2*d*x)^8 + 8*b^4*tan(1/2*d*x)^8 + 48*a^4*tan(1/
2*d*x)^7*tan(1/2*c) - 432*a^2*b^2*tan(1/2*d*x)^7*tan(1/2*c) + 96*b^4*tan(1/2*d*x)^7*tan(1/2*c) + 336*a^4*tan(1
/2*d*x)^6*tan(1/2*c)^2 - 2304*a^2*b^2*tan(1/2*d*x)^6*tan(1/2*c)^2 + 416*b^4*tan(1/2*d*x)^6*tan(1/2*c)^2 + 1200
*a^4*tan(1/2*d*x)^5*tan(1/2*c)^3 - 7344*a^2*b^2*tan(1/2*d*x)^5*tan(1/2*c)^3 + 1248*b^4*tan(1/2*d*x)^5*tan(1/2*
c)^3 + 2028*a^4*tan(1/2*d*x)^4*tan(1/2*c)^4 - 13392*a^2*b^2*tan(1/2*d*x)^4*tan(1/2*c)^4 + 2400*b^4*tan(1/2*d*x
)^4*tan(1/2*c)^4 + 1200*a^4*tan(1/2*d*x)^3*tan(1/2*c)^5 - 7344*a^2*b^2*tan(1/2*d*x)^3*tan(1/2*c)^5 + 1248*b^4*
tan(1/2*d*x)^3*tan(1/2*c)^5 + 336*a^4*tan(1/2*d*x)^2*tan(1/2*c)^6 - 2304*a^2*b^2*tan(1/2*d*x)^2*tan(1/2*c)^6 +
 416*b^4*tan(1/2*d*x)^2*tan(1/2*c)^6 + 48*a^4*tan(1/2*d*x)*tan(1/2*c)^7 - 432*a^2*b^2*tan(1/2*d*x)*tan(1/2*c)^
7 + 96*b^4*tan(1/2*d*x)*tan(1/2*c)^7 + 3*a^4*tan(1/2*c)^8 - 36*a^2*b^2*tan(1/2*c)^8 + 8*b^4*tan(1/2*c)^8 + 12*
a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + ta
n(1/2*c)^2 + 1))*tan(1/2*d*x)^6 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) +
2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)
^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6 - 12*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*
c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d
*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^6 + 18
*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/
2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + t
an(1/2*c)^2 + 1))*tan(1/2*d*x)^6 - 24*a^3*b*tan(1/2*d*x)^7 + 36*a*b^3*tan(1/2*d*x)^7 + 72*a^3*b*log(2*(tan(1/2
*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)
^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan
(1/2*d*x)^5*tan(1/2*c) - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/
2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1
/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c) - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1
/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/
2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*t
an(1/2*c) + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/
2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan
(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^5*tan(1/2*c) - 384*a^3*b*tan(1/2*d*x)^6*tan(1/2*c) + 432*a*b^3*t
an(1/2*d*x)^6*tan(1/2*c) - 2304*a^3*b*tan(1/2*d*x)^5*tan(1/2*c)^2 + 2256*a*b^3*tan(1/2*d*x)^5*tan(1/2*c)^2 - 4
56*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(
1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 +
 tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^3 + 684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x
)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) +
 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^3 + 456*a^3*b
*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2
*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^3 - 684*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(
1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan
(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c)^3 - 6000*a^3*b*tan(1/
2*d*x)^4*tan(1/2*c)^3 + 6120*a*b^3*tan(1/2*d*x)^4*tan(1/2*c)^3 - 6000*a^3*b*tan(1/2*d*x)^3*tan(1/2*c)^4 + 6120
*a*b^3*tan(1/2*d*x)^3*tan(1/2*c)^4 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c)
 + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d
*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^5 - 108*a*b^3*log(2*(tan(1/2*
d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^
2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(
1/2*d*x)*tan(1/2*c)^5 - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*
d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2
*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^5 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/
2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2
*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(
1/2*c)^5 - 2304*a^3*b*tan(1/2*d*x)^2*tan(1/2*c)^5 + 2256*a*b^3*tan(1/2*d*x)^2*tan(1/2*c)^5 + 12*a^3*b*log(2*(t
an(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(
1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1
))*tan(1/2*c)^6 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*t
an(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2
+ tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^6 - 12*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d
*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c)
 + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^6 + 18*a*b^3*log(2*(tan(1/
2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c
)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*ta
n(1/2*c)^6 - 384*a^3*b*tan(1/2*d*x)*tan(1/2*c)^6 + 432*a*b^3*tan(1/2*d*x)*tan(1/2*c)^6 - 24*a^3*b*tan(1/2*c)^7
 + 36*a*b^3*tan(1/2*c)^7 - 12*a^4*tan(1/2*d*x)^6 + 72*a^2*b^2*tan(1/2*d*x)^6 - 16*b^4*tan(1/2*d*x)^6 - 144*a^4
*tan(1/2*d*x)^5*tan(1/2*c) + 720*a^2*b^2*tan(1/2*d*x)^5*tan(1/2*c) - 96*b^4*tan(1/2*d*x)^5*tan(1/2*c) - 648*a^
4*tan(1/2*d*x)^4*tan(1/2*c)^2 + 3312*a^2*b^2*tan(1/2*d*x)^4*tan(1/2*c)^2 - 480*b^4*tan(1/2*d*x)^4*tan(1/2*c)^2
 - 1200*a^4*tan(1/2*d*x)^3*tan(1/2*c)^3 + 7344*a^2*b^2*tan(1/2*d*x)^3*tan(1/2*c)^3 - 1248*b^4*tan(1/2*d*x)^3*t
an(1/2*c)^3 - 648*a^4*tan(1/2*d*x)^2*tan(1/2*c)^4 + 3312*a^2*b^2*tan(1/2*d*x)^2*tan(1/2*c)^4 - 480*b^4*tan(1/2
*d*x)^2*tan(1/2*c)^4 - 144*a^4*tan(1/2*d*x)*tan(1/2*c)^5 + 720*a^2*b^2*tan(1/2*d*x)*tan(1/2*c)^5 - 96*b^4*tan(
1/2*d*x)*tan(1/2*c)^5 - 12*a^4*tan(1/2*c)^6 + 72*a^2*b^2*tan(1/2*c)^6 - 16*b^4*tan(1/2*c)^6 + 72*a^3*b*tan(1/2
*d*x)^5 - 60*a*b^3*tan(1/2*d*x)^5 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c)
+ 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*
x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c) - 108*a*b^3*log(2*(tan(1/2*d
*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2
 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1
/2*d*x)^3*tan(1/2*c) - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d
*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*
c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan(1/2*c) + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2
*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*
d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^3*tan
(1/2*c) + 720*a^3*b*tan(1/2*d*x)^4*tan(1/2*c) - 600*a*b^3*tan(1/2*d*x)^4*tan(1/2*c) + 312*a^3*b*log(2*(tan(1/2
*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)
^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan
(1/2*d*x)^2*tan(1/2*c)^2 - 468*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(
1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan
(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 - 312*a^3*b*log(2*(tan(1/2*d*x)^2*
tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*t
an(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x
)^2*tan(1/2*c)^2 + 468*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^
2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2*tan(1/2*c)^2 + 2304*a^3*b*tan(1/2*d*x)^3*tan(1/2*c)^2 -
 2256*a*b^3*tan(1/2*d*x)^3*tan(1/2*c)^2 + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1
/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(
1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^3 - 108*a*b^3*log(2*(tan
(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/
2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))
*tan(1/2*d*x)*tan(1/2*c)^3 - 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan
(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*ta
n(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c)^3 + 108*a*b^3*log(2*(tan(1/2*d*x)^2*t
an(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*ta
n(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)
*tan(1/2*c)^3 + 2304*a^3*b*tan(1/2*d*x)^2*tan(1/2*c)^3 - 2256*a*b^3*tan(1/2*d*x)^2*tan(1/2*c)^3 + 720*a^3*b*ta
n(1/2*d*x)*tan(1/2*c)^4 - 600*a*b^3*tan(1/2*d*x)*tan(1/2*c)^4 + 72*a^3*b*tan(1/2*c)^5 - 60*a*b^3*tan(1/2*c)^5
+ 18*a^4*tan(1/2*d*x)^4 - 72*a^2*b^2*tan(1/2*d*x)^4 + 144*a^4*tan(1/2*d*x)^3*tan(1/2*c) - 720*a^2*b^2*tan(1/2*
d*x)^3*tan(1/2*c) + 96*b^4*tan(1/2*d*x)^3*tan(1/2*c) + 336*a^4*tan(1/2*d*x)^2*tan(1/2*c)^2 - 2304*a^2*b^2*tan(
1/2*d*x)^2*tan(1/2*c)^2 + 416*b^4*tan(1/2*d*x)^2*tan(1/2*c)^2 + 144*a^4*tan(1/2*d*x)*tan(1/2*c)^3 - 720*a^2*b^
2*tan(1/2*d*x)*tan(1/2*c)^3 + 96*b^4*tan(1/2*d*x)*tan(1/2*c)^3 + 18*a^4*tan(1/2*c)^4 - 72*a^2*b^2*tan(1/2*c)^4
 - 12*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + t
an(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^
2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2 + 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2
*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/
2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2 + 12*a^3*b*log(2*(tan(1/2*d*x)^2*ta
n(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan
(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^
2 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 +
tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)
^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)^2 - 72*a^3*b*tan(1/2*d*x)^3 + 60*a*b^3*tan(1/2*d*x)^3 - 72*a^3*b*log(2*(t
an(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(
1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1
))*tan(1/2*d*x)*tan(1/2*c) + 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*ta
n(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*t
an(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c) + 72*a^3*b*log(2*(tan(1/2*d*x)^2*tan
(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(
1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*t
an(1/2*c) - 108*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/
2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan
(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*d*x)*tan(1/2*c) - 384*a^3*b*tan(1/2*d*x)^2*tan(1/2*c) + 432*a*b^3*tan
(1/2*d*x)^2*tan(1/2*c) - 12*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2
*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/
2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^2 + 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*ta
n(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan
(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^2 + 12*a^3*b*log(2*
(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + ta
n(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 +
 1))*tan(1/2*c)^2 - 18*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*tan(1/2*d*x)
*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^
2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1))*tan(1/2*c)^2 - 384*a^3*b*tan(1/2*d*x)*tan(1/2*c)^2 + 432*a*b^3*tan(1/2
*d*x)*tan(1/2*c)^2 - 72*a^3*b*tan(1/2*c)^3 + 60*a*b^3*tan(1/2*c)^3 - 12*a^4*tan(1/2*d*x)^2 + 72*a^2*b^2*tan(1/
2*d*x)^2 - 16*b^4*tan(1/2*d*x)^2 - 48*a^4*tan(1/2*d*x)*tan(1/2*c) + 432*a^2*b^2*tan(1/2*d*x)*tan(1/2*c) - 96*b
^4*tan(1/2*d*x)*tan(1/2*c) - 12*a^4*tan(1/2*c)^2 + 72*a^2*b^2*tan(1/2*c)^2 - 16*b^4*tan(1/2*c)^2 + 6*a^3*b*log
(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 +
 tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^
2 + 1)) - 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 + 2*tan(1/2*d*x)^2*tan(1/2*c) + 2*tan(1/2*d*x)*tan(1/2*c)
^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 - 2*tan(1/2*d*x) - 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2
*d*x)^2 + tan(1/2*c)^2 + 1)) - 6*a^3*b*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x)^2*tan(1/2*c) - 2*ta
n(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) + 1)/(tan(1/2*d*x)^2*t
an(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1)) + 9*a*b^3*log(2*(tan(1/2*d*x)^2*tan(1/2*c)^2 - 2*tan(1/2*d*x
)^2*tan(1/2*c) - 2*tan(1/2*d*x)*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 2*tan(1/2*d*x) + 2*tan(1/2*c) +
 1)/(tan(1/2*d*x)^2*tan(1/2*c)^2 + tan(1/2*d*x)^2 + tan(1/2*c)^2 + 1)) + 24*a^3*b*tan(1/2*d*x) - 36*a*b^3*tan(
1/2*d*x) + 24*a^3*b*tan(1/2*c) - 36*a*b^3*tan(1/2*c) + 3*a^4 - 36*a^2*b^2 + 8*b^4)/(d*tan(1/2*d*x)^8*tan(1/2*c
)^8 - 2*d*tan(1/2*d*x)^8*tan(1/2*c)^6 - 12*d*tan(1/2*d*x)^7*tan(1/2*c)^7 - 2*d*tan(1/2*d*x)^6*tan(1/2*c)^8 + 1
2*d*tan(1/2*d*x)^7*tan(1/2*c)^5 + 52*d*tan(1/2*d*x)^6*tan(1/2*c)^6 + 12*d*tan(1/2*d*x)^5*tan(1/2*c)^7 + 2*d*ta
n(1/2*d*x)^8*tan(1/2*c)^2 + 12*d*tan(1/2*d*x)^7*tan(1/2*c)^3 - 76*d*tan(1/2*d*x)^5*tan(1/2*c)^5 + 12*d*tan(1/2
*d*x)^3*tan(1/2*c)^7 + 2*d*tan(1/2*d*x)^2*tan(1/2*c)^8 - d*tan(1/2*d*x)^8 - 12*d*tan(1/2*d*x)^7*tan(1/2*c) - 5
2*d*tan(1/2*d*x)^6*tan(1/2*c)^2 - 76*d*tan(1/2*d*x)^5*tan(1/2*c)^3 - 76*d*tan(1/2*d*x)^3*tan(1/2*c)^5 - 52*d*t
an(1/2*d*x)^2*tan(1/2*c)^6 - 12*d*tan(1/2*d*x)*tan(1/2*c)^7 - d*tan(1/2*c)^8 + 2*d*tan(1/2*d*x)^6 + 12*d*tan(1
/2*d*x)^5*tan(1/2*c) - 76*d*tan(1/2*d*x)^3*tan(1/2*c)^3 + 12*d*tan(1/2*d*x)*tan(1/2*c)^5 + 2*d*tan(1/2*c)^6 +
12*d*tan(1/2*d*x)^3*tan(1/2*c) + 52*d*tan(1/2*d*x)^2*tan(1/2*c)^2 + 12*d*tan(1/2*d*x)*tan(1/2*c)^3 - 2*d*tan(1
/2*d*x)^2 - 12*d*tan(1/2*d*x)*tan(1/2*c) - 2*d*tan(1/2*c)^2 + d)

Mupad [B] (verification not implemented)

Time = 8.15 (sec) , antiderivative size = 268, normalized size of antiderivative = 1.49 \[ \int \sin (c+d x) (a+b \tan (c+d x))^4 \, dx=-\frac {{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^2\,\left (6\,a^4-48\,a^2\,b^2+\frac {32\,b^4}{3}\right )+2\,a^4\,{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^6+\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )\,\left (12\,a\,b^3-8\,a^3\,b\right )-2\,a^4-\frac {16\,b^4}{3}+24\,a^2\,b^2-{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^4\,\left (6\,a^4-24\,a^2\,b^2\right )-{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^7\,\left (12\,a\,b^3-8\,a^3\,b\right )-{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^3\,\left (20\,a\,b^3-24\,a^3\,b\right )+{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^5\,\left (20\,a\,b^3-24\,a^3\,b\right )}{d\,\left ({\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^8-2\,{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^6+2\,{\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )}^2-1\right )}-\frac {\mathrm {atanh}\left (\mathrm {tan}\left (\frac {c}{2}+\frac {d\,x}{2}\right )\right )\,\left (12\,a\,b^3-8\,a^3\,b\right )}{d} \]

[In]

int(sin(c + d*x)*(a + b*tan(c + d*x))^4,x)

[Out]

- (tan(c/2 + (d*x)/2)^2*(6*a^4 + (32*b^4)/3 - 48*a^2*b^2) + 2*a^4*tan(c/2 + (d*x)/2)^6 + tan(c/2 + (d*x)/2)*(1
2*a*b^3 - 8*a^3*b) - 2*a^4 - (16*b^4)/3 + 24*a^2*b^2 - tan(c/2 + (d*x)/2)^4*(6*a^4 - 24*a^2*b^2) - tan(c/2 + (
d*x)/2)^7*(12*a*b^3 - 8*a^3*b) - tan(c/2 + (d*x)/2)^3*(20*a*b^3 - 24*a^3*b) + tan(c/2 + (d*x)/2)^5*(20*a*b^3 -
 24*a^3*b))/(d*(2*tan(c/2 + (d*x)/2)^2 - 2*tan(c/2 + (d*x)/2)^6 + tan(c/2 + (d*x)/2)^8 - 1)) - (atanh(tan(c/2
+ (d*x)/2))*(12*a*b^3 - 8*a^3*b))/d